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ABSTRACT 

Development of a new and efficient sub-band coding of digital 
images is reported This includes 1) the design of short kernel 
symmetric analysis/synthesis filter banks with perfect reconstruction 
capability; and 2) entropy coding of the sub-bands by using 
arithmetic coding in combination with DPCM coding of the lowest 
band and PCM coding of the higher bands. 



1. INTRODUCTION 

Subband coding [I J, in general, refers to the class of coding 
techniques where by parallel application of a set of filters (i.e., 
analysis filters) the input signal is decomposed into several 

' narrowbands. " The resulting bands are then decimated and coded 
separately for the purpose of transmission. For reconstruction, 
decimated signals are decoded, interpolated, and filtered (i.e; 

^synthesis filtering) before being added to reproduce the original 
signal.: . '"This approach, in general, demands the design of 

• sophisticated band-pass filters to minimize the perception of aliasing 
effects, jn a much simpler approach, Esteban et al. [2] introduced the 
concept of .Quadrature Mirror Filters (QMFs) which, in the absence 
of channel .and quantization noise, permits an alias-free and near 
perfect reconstruction of the input signal. Theoretical extension of 
dnc-dimensional (1-D) QMF filtering to multi-dimensional cases has 

';bcen treated by Vetterli [3J. Moreover, application of two- 
; dimensional (2-D) QMF concept to the sub-band coding of digital 
images has recently received considerable attention and it has been 
shown to be highly effective for image compression [4J-[8] 

A disadvantage with the classical QMF approach, as described in [2J. 
- is that it docs not permit reconstruction to be exact, although, 

• amplitude distortion can be made very small by using long tap filters. 
; For video arid digital image applications use of such long tap filters, 

• while not providing any significant coding gain, may increase the 
hardware complexity. 

.jln this paper; analysis/synthesis filter banks arc realized by using 
'symmetric s'hort kernel filters with perfect reconstruction capability. 
Furthermore, an adaptive sub-band coding scheme which uses 
arithmetic coding in combination with the DPCM coding of the 
lowest band and PCM coding of upper band signals is also presented 
In Section II principle of filtering operation is discussed Our 
proposed coding scheme is then described in Section 3. Simulation 
results are covered in Section 4. Summary and conclusions appear in 
Section 5. 



1 Filtering: Exact Reconstruction Filter banks 
Figure 1 depicts the principle elements of a two-channel 
analysis/synthesis filter bank system. As shown, the input signal *fn), 
is filtered into two bands, which are then decimated by a factor of 2. 
At the synthesis side, the signals ™(«) and are interpolated (Le4 
by inserting zeros between successive samples), filtered and added 
together to reproduce the origir^ signaL Output Y{x) in this figure 
can be written as: 



y{x) - g 0 [z)y 0 (*) + Gi(OriM 
where 

y 0 (*) - f [*(>(*)*"(*) + »6<-.*)f.i-*)1 : 

rid) - jWib)X(t) +: //i(-i)jf(-x)l > 
Substitution of Eqs. (22H2.3) into (2.1) yields 
y(x)«i/ 2 fG 0 (x)// D (*) + thCO"i(i>l*^V 

+ V2lC Q iz)H 0 i-i)+ G l (x)tf l (-*)]f(--xj 



(2.1) 

(2.2) 
(23) 

(2-4) 



From Eq. (Z4) it is easily seen : that aliasing is removed if the 
synthesis filters are defined as 



C 0 (i) - Hxi-j) 
Gi(x) = - Hq(-i) 



(2.5.1) 
(2.5.2) 



For this class of synthesis .filters, perfect reconstruction of the input 
signal requires that I;;- 



A(x) = />(i)- - ax- 
where 

Pit) - HoWHii-x). 



(2.6) 



(2.7) 



is called the product filter, and m is . the delay introduced by the 
combination of analysis/synthesis filter banks 

11 Design of Short Kernel Filter Bank 

The classical approachlto the design 'of analysis/synthesis filter banks 
which removes aliasing, and'- approximates Eq. (Z6) is based on the 
Quadrature Mirror Filter (QMF) concept [2j. This approach makes 
use of even, symmetrical FIR -filters. In addition, QMFs are defined 
to be mirror image' of one another i.e^ «i(x) - H 0 {-x). For this class 
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of analysis/synthesis filter banks it is not possible to obtain exact 
reconstruction if, FIR filters greater than length two are used As an 
alternative, Smith and Barnwell [9] have proposed a design technique 
which consists of designing the product filter P(z), and then factoring 
the result into separate analysis synthesis filter banks (e.g.; minimum 
and maximum phase components). This technique, in principle, does 
not lead to linear phase analysis/synthesis filters, and it may not be 
suitable for image coding applications where the use of symmetric 
filters is more desirable. 

In this paper we address the design of symmetric short tap filters. 
Our design procedure is based on the factorization of product filter 
j»U) into linear phase components as suggested by Vettcrli 110} 
Referring to Eqs. (Z6)-(2.7), a number of observations can be made. 
First, from Eq. (27) note that if both Hq{x) and Hi(-i) arc- linear 
phase filters, then their product P(z) also becomes a linear phase 
filter. Second, Eq. (2.6) can be interpreted as a constraint on the 
coefficients of P(z). It then follows that product filters with an odd 
number of coefficients can only satisfy (26). 

Our design strategy is to choose a half-band product filter P(z) as 
shown below 

P[x) = co + a 2 j " 2 +' . . . + 02,-2* -2p+2 + * 

+ a 2p< *' 7p+2 + ....... + ooz . . \ . (3.1) 

and optimize the coefTicients based on a series of constraints on the 
location of the zeros. Note that for a P(z) of length 4p-l, there are 
only p coefTicients to be optimized and the result will yield two 
filters of length 2p. Among the design constraints one can impose is 
to have filters with integer power of two coefficients. 

Case p =1; 

This is the simplest filter bank and the form of P(z) is 

i»(z) - oo + ^+aor" 2 - (3.2) 

The only solution that is a low pass filter is: . 

P(x) = (1/2 + 1/22 7*) (1/2 + 1/21-1). (33) 

This is a unique linear phase filter bank of length- two and it is the 
only QMF solution with exact reconstruction as well (see [9]). Such 
filters can only be useful if coding accuracy is high, otherwise the 
type of artifacts introduced by quantization noise makes the resulting 
picture visually unpleasant 

Case p = 2; 

In this case the form of P(z) is . 

P(i) = a 0 + « 2 i~ 2 . + *~ 3 + + a Q*~* ■ (3.4) 

One degree of freedom can be easily suppressed by imposing a 
(double) zero at z - - 1 . P(z) can then be factorized as 

+ (1/1 + 2fl 0 ) i" 2 - 2a 0 z" 3 + a 0 x "«). (3.5) 

Note that useful filters exist only for values of a 0 in the interval [- 
1/8, 0J. For example, by setting * a 0 ~-vi6 we get an additional 
double zero at -1. In this case, the expression for P(z) becomes 

/»(*) = 1/16 (l + r-l) 3 (-1 + 3X-1 +5r -2- . (3.6) 

This factorization yields two pairs of filters for exact reconstruction 



Solution 1 (equal length): 

Wo(0 = + si -1 + 3j~ 2 + i"" 3 ) 

Ni(-*) « 1/4 (-1 + si" 1 + ai" 2 - / " 3 ) (3.12) 

J/fl(x) - V4 (-1 + Si ' l + 3J ' 2 " I ~ 3 ) ( 3 ' 7 '^ 

H,(-z) - 1/4(1 + 3x" 1 + 3*~ 2 + x' 3 ) (3- 7 - 4 ) 

Solution 2 (unequal- length): 

H 0 {z) = 1/8 (-1 + 22-1 + 6,-2 + ^-3^-4) (3.8.1) 

' o 1/2 (1 + 21-! + m<) (3.8^) 



By comparing the filters in solution 1 and solution 2, it appears that 
the interpolation process to reconstruct the baseband signal in 
method 2 is what is commonly known as "linear interpolation*. This 
interpolation scheme is known to produce visually pleasant and 
smooth outputs. From the point of view of computational 
complexity both filter banks can be implemented by a limited 
number of shifts and adds, and thus are computationally more 
efficient than their traditional QMF counterparts. 

3. CODING * V 

For channel coding of sub-bands, the basic .approach of DPCM. 
coding of lowest band signal, PCM coding of upper band signals, and 
run-length coding of the address of non-zero PCM values has been 
reported by [Gharavi, and Tabatabai I986J. In order to achieve a 
higher compression rate, we have extended the above method by 
using arithmetic coding [11] in combination with DPCM coding of 
lower band and PCM coding of higher band~ signals. Arithmetic 
coding method used in. our simulations is referred to as "Q-Coder" 
and has been described in detail by [Mitchell, and Pennebaker 1987]. 
The Q-coder has several important . properties; it avoids 
multiplication, and it does not require any predetermined . symbol 
probabilities. The principle characteristics of our coding method is 
described in the following section. " 

3.1 Coding of the Lowest Band Signal 1 

Similar "to the approach taken in [13], we have combined adaptive 
DPCM and arithmetic coding. ' Briefly/in. this technique, a Markov- 
like state is determined from a measure of local image activity, and 
the Q-coder is used to encode the magnitude of the difference data \ 
by using . probabilities , conditioned on. the state. Separately, 
conditional entropy coding is also applied to encode the sign of error 
data. Here, for our purposes, we have used a different set of 
quantization levels and threshold values. This is because the quality 
of the reconstructed, picture depends, to large extent, on the fidelity 
with which the lowest band signal is coded. Tables I and II, shown 
below, illustrate the quantizer levels and their corresponding 
threshold values. 

In Table !, for example, all differences from 8 to 10 are mapped to 9. 

3.2 Coding of the Higher Bands 

Our basic strategy , to encode higher band signals is based on using a 
symmetric PCM quantizer with the following characteristics [5f 1) it 
contains a dead- zone, d; 2) an active range where signals falling in 
this range are uniformly coded into L levels; and 3) a saturation value 
Y, for signals whose magnitude are above a given threshold l 

By using the above quantizer a binary sequence suitable for the input 
to the Q-coder is generated. That is, first depending on whether. the 
magnitude of the signal is below or above "d" a M 0" or "1" is sent 
Note that such a binary sequence, in general, will correspond to the 
contour information of the higher band signals. For cases when a "1" 
• is sent a code word for the sign and magnitude of the non-zero PCM 
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TABLES I-II: QUANTIZATION AND THRESHOLD VALUES s SUMMARY AND CONCLUSIONS * 



Table! 




Tfthle 


II 


Threshold 


Level 




Threshold 


Level 


5 


3 




11 


6 


8 


6 




21 


15 


11 


9 




31 


25 


15 


12 




41 


35 


20 


17 




52 


45 


26 


22 




72 


60 


. 34 


29 




95 


83 


44 


38 




180 


135 


57 


49 




256 


255 


74 


64 






95 


83 






121 


106 






175 


135 






256 


255 







In this paper a simple, and efficient method of sub-band coding of 
digital images has been reported. First, a technique for designing 
symmetric short tap filters has been presented, and it has been shown 
that such filters can be easily; implemented by using simple 
arithmetic operations (eg.; addition and multiplication). By applying 
the above filters the input image has been decomposed into four 
bands, which have then been coded by using arithmetic coding in 
combination with DPCM coding of the lowest band and PCM coding 
of higher bands. 

Our simulation results, have demohstrated that by using the method 
mentioned above good quality, pictures can be obtained in the ranee 
of 0.7 to 0.8 bits/pel. • 
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Table III -- Bit Rate and SNR for LENA and PEPPERS 



Picture 


Bit Rate 


SNR 
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.7 


34.56 


PEPPERS 


.77 


33.83 
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